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The earlier-developed master equation approach and kinetic cluster methods are applied to study 
kinetics of Llo-type orderings in alloys, including the formation of twinned structures characteristic 
of cubic-tetragonal-type phase transitions. A microscopical model of interatomic deformational 
interactions is suggested which generalizes a similar model of Khachaturyan for dilute alloys to 
the physically interesting case of concentrated alloys. The model is used to simulate Al^Llo 
transformations after a quench of an alloy from the disordered Al phase to the single-phase Llo 
state for a number of alloy models with different chemical interactions, temperatures, concentrations, 
and tetragonal distortions. We find a number of peculiar features in both transient microstructures 
and transformation kinetics, many of them agreeng well with experimental data. The simulations 
also demonstrate a phenomenon of an interaction-dependent alignment of antiphase boundaries in 
nearly-equilibrium twinned bands which seems to be observed in some experiments. 

I. INTRODUCTION 

Studies of microstructural evolution under alloy phase transformations from the disordered FCC phase (Al phase) 
to the CuAu Ttype ordered tetragonal phase (Llo phase) attract interest from both fundamental and applied points 
of view. A characteristic feature of such transitions is the formation in the ordered phase of peculiar 'polytwinned' 
structures consisting of arrays of ordered bands separated by the antiphase boundaries (APBs) lying in the (llO)-type 
planes, while the tetragonal axes,of antiphase-ordered domains (APDs) in the adjacent bands have 'twin-related' 
(100) and (OlO)-type orientationsETLl Transformation Al^Llo includes a number of intermediate stages, including 
the 'tweed' stage discussed below. These transformations are inherent, in particular, to many alloy systems with out- 
standing magnetic characteristics, such as Co-Pt, Fe-Pt, Fe-Pd and similar alloys, and studies of their microstructural 
features, for example, properties and evolution of APBs, are interesting for applications of these systems in various 
magnetic devices for which the structure and the distribution of APBs can be very importantErti. n ._. 

The physical reason for the formation of twinned structures was discussed by a number of authorscltlJ, and it is 
explained by the elimination of the volume-dependent part of elastic energy for such structures. However, theoretical 
treatments of the kinetics of Al— >Llo transformation seem to be rather scarce as yet. Khachaturyan and coworkerstj 
discussed kinetics of tweed and twin formation using a 2D model in a square lattice with a number of simplifying 
approximations: a mean-field-type kinetic equation; a phenomenological description of interaction between elastic 
strains and local order parameters; an isotropic elasticity; an unrealistic interatomic interaction model (with the 
nearest- neighbour interaction being by an order of magnitude weaker than more distant interactions), etc. In spite 
of all these assumptions, seme features of evolution found by Khachaturyan and coworkerscj agree qualitatively with 
experimental observationsBa. It may illustrate a low sensitivity of these features to the real structure and interactions 
in an alloy. However, such an oversimplified approach is evidently insufficient to study the details of evolution and their 
dependence on the characteristics of an alloy, such as the type of interatomic interaction, concentration, temperature, 
etc, which seems to be most interesting for both applications and physical studies of the problem. 

In this work we investigate kineticSj-of the Al— >Llo transition using the microscopical master equation approach 
and the kinetic cluster field methodE3E3. Earlier this method was used to study Al-^Ll2-typc transformationsE3 as 
well as early stages of the Al^Llo transition when the-,deformational interaction Hd due to the tetragonal distortion 
of the Llo phase is still insignificant for the evolutionlla. Here we consider all stages of this transition, including the 
tweed and twin stages when the interaction Hd becomes important. To this end w first derive a microscopical model 
for Hd which generalizes the analogous model of Khachaturyan for dilute alloysli3 to the physically interesting case 
of concentrated alloys. Then we employ the kinetic cluster field method to simulate Al— >Llo transformation in the 
presence of deformational interaction Hd for a number of alloy models with both short-range and extended-range 
chemical interactions at different temperatures, concentrations and tetragonal deformations. The simulations repeal 
a number of interesting microstructural features, many of them agreeing well with experimental observationsaQ. 
We observe, in particular, a peculiar phenomenon of an interactiq&rdependent alignment of orientations of APBs 
within twin bands which was earlier discussed phenomenologicallyliil. The simulations also show that the type of 
microstructural evolution strongly depends on the interaction type as well as on the concentration c and temperature 
T . In particular, drastic, phase-transiton-like changes in morphology of APBs within twin bands can occur under 
variation of c or T in the short-range-interaction systems. 
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The paper is organized as follows. In section 2 we derive a microscopical expression for the deformational interaction 
Hd in concentrated alLws. In section 3 we describe our methods of simulation of Al— >Llo transition which are similar 
to those used earlieiOlla. In section 4 we investigate the transformation kinetics for the alloy systems with an extended 
or intermediate interaction range, and in section 5, that for the short-range-interaction systems. Our main conclusions 
are summarized in section 6. 



II. MODEL FOR DEFORMATIONAL INTERACTION IN CONCENTRATED ALLOYS 

We consider a binary substitutional alloy A c Bi_ c . Various distributions of atoms over lattice sites i are described 
by the sets of occupation numbers {rii} where the operator Uj = n^i is unity when the site i is occupied by atom A 
and zero otherwise. The effective Hamiltonian H e f{ describing the energy of these distributions has the form 

H e ff = 2J v lj n l n j + ii,n .us, + . . . (1) 

i>3 i>j>k 

where are effective interactions. 

The interactions include the 'chemical' contributions v% j which describe the energy changes under the sub- 
stitution of some atoms A by atoms B in the rigid lattice, and the 'deformational' interactions vf j related to the 
difference in the lattice deformation under such a substitution. A microscopical model for v d in dilute alloys was 
suggested by Khachaturyar£3. The deformational interaction in concentrated alloys can lead to some new effects that 
are absent in the dilute alloys, in particular, to the lattice symmetry changes under phase transformations, such as 
the tetragonal distortion under Ll ordering. Earlier these effects were treated only phenomenologicallytJ. Below we 
describe a microscopical model for calculations of v d which generalizes the Khachaturyan's approachEJ to the case of 
concentrated alloys. 

Let us denote the position of site k in the disordered 'averaged' crystal as r^. Because of the randomness of a 
real disordered or partially ordered alloy the actual atomic position (averaged over thermal vibrations) is not but 
rfe +Ufe where uj. is the 'static displacement'. Supposing this displacement to be small we can expand the 'adiabatic' 
(averaged over rapid phonon motion) alloy energy H — H{rii, u^.} to second order in uj,: 

H = H c {rii} - U ak F ak + - u akUf3lA a k,/3l (2) 
k k,l 

where a and are Cartesian indices and the summation over repeated Greek indices is implied here and below. The 
term if c {n,} in (^) describes interactions in the undistorted average crystal lattice, i. e. chemical interactions v\ . 
mentioned above; F ak can be called the generalized Kanzaki force; and A a k } pi is the force constant matrix. Both 
quantities F ak and Aak,pi are certain functions of occupation numbers and their evaluation needs some further 
approximations. 

Below we consider ordering phase transitions at the fixed-mean concentration c. Changes of elastic constants 
and phonon spectra under such transitions are usually smalltll. Therefore, the force constant matrix A a k,pi can be 
reasonably well approximated with the simple 'average crystal' approximation: A a k t pi{rii} — ► A a k t pi{c} = A a j~,pi- To 
approximate the Kanzaki force F ak we first formally write it as a series in the occupation numbers nf. 

F ak {n t } = F« k + ]T F^m + ]T F^ ni n 3 +... (3) 

i i>j 

Equilibrium values of displacements = u^{rii} at the given distribution {n^} are determined by the minimization 
of energy (||) over uj-, and the constant F^ k in ([|) affects only the reference point u|{0} in the function u|{nj}. 
This constant can be determined, for example, from the condition of vanishing of mean static displacements in the 
averaged crystal at some c = Co, which implies the relation: (F ak {rii}) c=Ca = where the symbol (. . .) means the 
statistical averaging over an alloy. The constants i^ fc are insignificant for what follows, and below they are omitted 
to simplify formulas. 

In writing an explicit expression for the contribution Hk (to be called for brevity the 'Kanzaki term') of the 
occupation-dependent Kanzaki forces in energy (p|) one should consider that due to the translation invariance it can 
include only differences of displacements (life — u^), (u^ — Uj), etc. Therefore, this term should have the form 



(4) 



where f^* 1 •••*"> = f*= are some parameters describing interaction of lattice deformations with site occupations. 

Representation (|J) for Hk as a sum of contributions of m-site 'clusters' proportional to products . . . Ui m is 
analogous to similar cluster expansions for the 'chemical' Hamiltonian H c {rii} in (||). These expansions have been 
widely discussed, in particular, in connection with first-principle calculations of chemical interactions u| i = w c ' m , 
see e. g.lla. The calculations have shown that the values of m-site interactions v c ' m in most alloys rapidly decrease 
with an increase of to, and the pairwise interaction v c ' 2 is usually dominant. It is natural to expect that a similar rapid 
convergence is also typical for the expansion (Q). Therefore, below we omit many-site interactions with to > 2 in 
Eq. (Q). At the same time, in estimates of parameters for real alloys below we combine some model assumptions 
about with using of available experimental data about the variations of lattice deformations with concentration 
and orderings, and such estimates may also implicitly include the contributions of many-site interactions f^. 

For what follows it is convenient to proceed from functions Uj, = u(r/ £ ), rii — n(rj), fj'' = fi(rfc — r^), i 2 ' lJ = 
ki^k — r i> r j — r i) an d Aak^pi — A a p(rk — Yi) in Eqs. (|J) and (ji|) to their Fourier components in the average crystal 
lattice. Then the energy (g) takes the form: 



H = H c { ni } + 1 u -k ( "kfik + E CT k f 2k ) + E 

k V R / k 



u^Aful (5) 



Here N is the total number of crystal cells, the summation over k goes within the Brillouin zone of the averaged 
crystal, and we use the following notation: 

u k = ]T u(r)e-* kr ; n k = ^ n(r)e~ ikr ; o* = ^ n(r)n(r + R)e" ikr ; 

r r r 

f lk = 5Xr)(l-e- ikr ); f* = ]T f 2 (r, R)(l - e^)- ^^I^(r) e - b (6) 

r r r 

If one adopts a commonly used model of 'central' Kanzaki forces in which forces f-f" 1 and f 2 ' < ' 5 in (|]) are supposed to 
be proportional to the vector = (r^ — r^), the vector functions fik and f2k m (m can be expressed via two scalar 
functions, ifi and if2'- 

f lk = ^r^( r )(l-e- k '), f^^r^RKl-^). (7) 

r r 

The functions ipi and ip 2 in (^) determine the dependence of equilibrium lattice parameters on concentration or 
ordering. To show it we first note that the homogeneous deformation u a p is described by Fourier-components Uk 
with small k — * 0, while functions fik and f2k in Eqs. (||) and (|J) at small k are linear in k. Thus the contribution of 
homogeneous deformations to the Kanzaki term in (^) is proportional to Fourier-components u^P of the elastic strain 
u a /3 = (du a /dx/3 + dup/dx a )/2 at k — > and, according to first equation (0), these components are related to u a p as 



a/3 



| k ^o = i{kpv£ + k a u£)\ k ^ Q = Nu a p. (8) 



At small k the force constant matrix A k in is bilinear in k, and the last term of (||) corresponds to the standard 
expression for the elastic energy bilinear in u a p and linear in the elastic constants c a p^s, see e.g.tl Therefore, the 
total contribution of terms with the homogeneous elastic strain u a p to energy (|^) (to be called 'the elastic strain 
energy' E e i) can be written as 



Eei = -Uap Mi' 3 n + E A^af^j + i Nil, c afil& u a 



pUyS- (9) 

Here f2 is the volume per atom in the average crystal; quantities A"^ and ^r are expressed via functions ipi and (p 2 
in as: 

= y^ x a x P ( Pi( r )> A 2R = ^2x a x ip 2 (r,Il), (10) 

r r 

where x a is the Cartesian component of vector r = (xi, x 2 , x^); and no or cr R is the Fourier component, rik or <r R , at 
k = 0. According to Eq. (||), the operator no or <r R is the sum of a macroscopically large number N of similar terms. 
Thus within the statistical accuracy each of these operators can be substituted by its average value: 
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n = N(n(r)) = Nc; <j r = N(n(r) n(r + R)). 



(11) 



The last average in (|ll|) can be expressed via mean occupations of sites and their correlators. In an ordered alloy 
there exist several non-equivalent sublattices s with the lattice vectors r s and mean occupations c s — (n(r s )), and so 
the last average in (Thy) includes averaging over all sublattices s: 



(n(r) n(r + R)) = ^ v s (c s c sR + K sR ) . 



(12) 



Here c s r is the mean occupation (n(r)) for r = r s + R; v s — N s /N is the relative number of sites in the sublattice s; 
and K sR is the correlator of occupations of sites located at r = r s and at r = r s + R: 



K sR = ([n(r g ) - c s ] [n(r a + R) - c sR ]). 



(13) 



In a disordered alloy all sites are equivalent, thus c s = c s r = c; u s = 1; and both index s and the summation over s 
in (12) are omitted. 

Using Eqs. ( |ll|) and (|lj) one can rewrite the elastic strain energy (||) as 



A°f C + ^ X! Vs ( Cs CsR + ^ sR ) A 2l 




2R 



R s 



-A^O C a f3 7 SU a f3U 7 S- 



(14) 



The correlator i4T S R in Eq. (14) can be calculated using that or another method of statistical theory. However, for 



most alloy systems of practical interest, in particular, at c and T values not close to the thermodynamic instabilty 
points T s , the correlators K sR are small and can be neglected. Then equation dl3) is simplified: 



E e i — —u<xp 



NAfc + J2<r)c(r + R)A^ 



r.R 



(15) 



Equilibrium values of u a p in the absence of applied stress are determined by the minimization of energy E e i with 
respect to u a p which gives: 



fi Ca^sujs = A°f c + — c(r)c(r + K)A. 



a/3 
2R- 



(16) 



r,R 



Eq. ( |l6| ) enables one to express the equilibrium strain u a p via the concentration, order parameters, and the interaction 
parameters A"^ and A^ R , and it can also be used to estimate these interaction parameters from experimental data 
on u a p(c,T). 

Let us consider Eqs. ( fl5| ) and ( fl£| ) in particular cases. For a disordered phase with c(r) = c, Eq. ( |l6| ) takes the 
form 



fi c af 3 7 sujs = A" 13 c + A% p c 



(17) 



where A^ — /Cr^2R- H ^ e disordered phase has a cubic symmetry (as for the FCC or BCC alloys), quantities 
A" and A£ are proportional to the Kronecker symbol S a /3, and Eq. ( |l7| ) determines the concentrational dilatation 
u(c) = u aa (c) - u aa (0): 



u{c) = (A lC + A 2 c 2 )/nB. 



(18) 



Here B = (en + 2ci2)/3 is the bulk modulus; Cij are the elastic constants in Voigt's notation; and coefficients Ai and 
A2 are expressed via functions tpi and 922 in (0), (|l0|): 



^=2^(^/3; A 2 =^> 2 (r,R)r 2 /3. 



(19) 



R 



The linear in c term in (|18|) corresponds to the Vegard law while the term with A2 describes the non-linear deviations 
from this law. Such deviations were observed for many alloys, and these data can be used to estimate A2 values, but 
in these estimates one should also take into consideration a possible concentration dependence of the bulk modulus 
B. 
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For the ordered phase, the mean occupation- c(r) can be written as a superposition of concentration waves corre- 
sponding to certain superstructure vectors k p E£l: 

c(r) = c+^[)jp exp(ik p r) + 77* exp(-ik p r)] , (20) 
p 

and amplitudes ??„ can be considered as order parameters. After the substitution of expressions ( BCj ) for c(r) and 
c(r + R) in Eq. ( |l5| ) the linear in ?7 P terms vanish due to the crystal symmetry, and the first term of (|l^) becomes 
the sum of the ordering-independent term and the term bilinear in order parameters: 

Eel = -Nu af j ^A" P C + Alf C 2 + ^ QaPpplVpl 2 ^ + ^NfLCaPjSUapUjS- (21) 

Here quantities q a f3 PP have a different form in the cases (a) when the superstructure vector k p is half of some reciprocal 
lattice vector g and thus both the order parameter r\ p and all factors exp(ik p r) in ( |20| ) are real, and (6) when k p ^ g/2: 

(a) k p = g/2 : q af3pp = ^ x a xp<p 2 {r, R) exp («k p r); (22) 

r : R 



(6) k p ^g/2: q a(ipp = ]-^x a xp(p 2 {r, R)cos(k p r). (23) 



2 

r,R 



The coefficients q a ji PP in ( pl| ) (to be called the 'striction' coefficients, in an analogy with the terminology used in the 
ferroelectijicity or magnetism theory) are commonly used in phcnomcnological theories of lattice distortions under 
orderingsail3. Eqs. (p2), (23) and ([LOh provide the microscopic expression for these coefficients via the function ip 2 
describing non-pairwise Kanzaki forces in Eqs. (^)-(0). 

Let us apply Eqs_/0)-(22) to the case of Llo or LI2 ordering in FCC alloys which are described by three real 
order parameters ry Q Ed'H Eqs. (|2C]) here take the form 

c(r) = c + r]i exp(ikir) + 772 exp(ik 2 r) + 773 exp(ik 3 r), (24) 

where k Q = g Q /2 is the superstructure vector corresponding to 770,: 

k x = [100]27r/a, k 2 = [010]2vr/a, k 3 = [001]27r/a. (25) 

In the cubic LI2 structure one has: |?7i| = \if 2 \ = ViVzVs > 0: an d f° ur types of ordered domains are possible. 
In the Llo-ordered structure with the tetragonal axis a a single parameter 770, is present which is either positive or 
negative, and so six types of ordered domains are possible. 

The striction coefficients for Llo or LI2 ordering are determined by Eq. (^2|). Due to the cubic symmetry of the 
'average' FCC crystal, there are only two different striction coefficients, qnn and qn 22 (and those obtained from 
them by the cubic symmetry operations), which for brevity will be denoted as qn and q± 2 , respectively: 

911 = ^Xi<£ 2 (r,R)exp(ikiR); q 12 = x 2 (p 2 (r, R) exp(ik 2 R). (26) 

r,R r,R 

Variation of elastic constants c a p 7 s with ordering is usually smallEll, and for simplicity it will be neglected. Then 
minimizing energy (|2l]) with respect to u a/ 3 we obtain the expressions for lattice deformations induced by ordering 

u = q+(ri 2 i+vt+vl)/nc + , e = q-[r&- (rj 2 + £ = q_(r% - r&)/Qc-. (27) 

Here u = u\\ +u 22 +U33 describes the volume change; e = u\x — (u 22 +U3a)/2 is the tetragonal distortion; £ = U22 — U33 
is the shear deformation; and q± or c± are linear combinations of striction or elastic constants: 

q-=qn~qi 2 ; c_ = c u -ci 2 ; q+ = qu + 2q 12 ; c + = cn + 2c 12 . (28) 

For the LI2 ordering, values | -771 1 = \r] 2 \ — \rj3\ — rj are the same, so just the volume striction u = 3q+rj 2 /£lc+ is 
present, while in the Lip-ordered domain with ?? 2 = 773 = one has both the volume and the tetragonal striction: 
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u = q + r] 2 /Qc + , e = q-ril/Qc-. (29) 

Therefore, using experimental data about the lattice distortions and order parameters under LI2 and Llo orderings 
one can estimate the striction coefficients qn and q± 2 and thus the non-pairwise Kanzaki interaction (p 2 in Eqs. (|2q). 

Below we suppose for simplicity the interaction (p 2 (r, R) to be short-ranged, i. e. significant only when each of 
three relative distances r, R and |r — R| does not exceed the nearest-neighbour distance p = a/y/2. Then this function 
can be written as 

Lp 2 (r, R) = S r ,pS RtP (ifia 5|r-R|,0 + 'fib S\ r -R\.p) (30) 

where <5 r . p is the Kronecker symbol equal to unity when r — p and zero otherwise while <p a and ifb are the interaction 
parameters. The assumption (|30| ) is analogous to that used by Khachaturyar£3 for the pairwise Kanzaki interaction 
Mr) in ©: 

M r ) = Mr,p (31) 

where the constant ip± is estimate*Lfrom experimental data on concentrational dilatation. First-principle estimates of 
lattice distortions in dilute alloystia seem to imply that the assumption ( |3l|) yields the correct order of magnitude of 
<£>i(r). Therefore, the analogous assumption (|30| ) for ^(r, R) can be reasonable, too. 

Substituting Eq. ( |30| ) into ( p6| ) we obtain the explicit expression for coefficients qik via parameters tp a and ipt in 
©: 

qu = -2a 2 (p a ; q 12 = -4a 2 ip b . (32) 
The coefficient in Eqs. (|1) and (0) for model ( |3l| ) has the foriiS 

3 

fik = 4:ipii ^ a a sin (ka a ) \_. cos O&p) (33) 

a=l P^a 

where a a is e a a/2 and e a is the unit vector along the main crystal axis a. The function in (||), (||) for model ( |30| ) 
is the sum of two terms: 

f 2 ^=fa R k+f 6 k- (34) 

Here, is ip a SR, p R (l ~ e~ jkR ), while the function for R equal to R a + R^ (where R Q is & a or (— a a ), R^ is 
a./3 or (— a^), and ^ a) can be written as 

= 2(p b [R + (ia 7 sin ka' - R Q cos ka') e ^ kR " + (ia' sin ka' - R /3 cos ka') e - 4kR f*] (35) 

where a' is [e Q e^]a/2. 

Relations (Q), j3C|)-(|35|) together with (|l8| ) and (|32] ) provide a simplified model for the Kanzaki term Hk in Eqs. (H) 
and (|5|). This model will be used below in simulations of A1^L1 transitions. To get an idea about the actual scale 
of parameters of this model, let us estimate quantities qik, <p a and ip^ in Eqs. (|3^) and (fj2^ fpr the alloys Co-Pt for 
which detailed data about the lattice distortion under LI2 and Llo orderings are availableE3£lL The volume change u 
under both the LI2 ordering in CoPt3 and Llo ordering in CoPt appears to be very smalEl'EIffil: u < 10 -3 . According 
to Eqs. ( p7|) and (|29| ) it implies the relationj-|<7i 2 ~ qu/2. The value q_ = (qu — qi 2 ) for CoPt can be estimated 
from second equation (E9I) using data of RefEl for 7/1 and e at T = 0.84 T c : r?i ~ 0.4; e ~ -0.04 (with the thermal 

expansion effect subtracted); and for the atomic volume: Vl — f2(ZL+ 0) ~ 13. 8A . Using also for the elastic constant 
c_ = (en — C12) its value for the FCC platinum, c_ ~ 0.97 MbarE 2 ], we obtain: g_ ~ 2.6 • 10 4 K. Combining it with 
the above-mentioned relation qi 2 ~ <7n/2 and using Eq. ( p7| ) we find: ip a ~ 2 • 10 4 K/a 2 , and </3fe ~ 5 • 10 3 K/a 2 . Let 
us also note that the ordering-induccd clastic energy per atom e°\ d in the CoPt alloy is small: e°J d ~ f2c~£ 2 /6 ~ 
30 K, which is much less than the Llo ordering temperature T c ~ 1100 K. 

The equilibrium values of displacements u£ = u^.(ni) are found by the minimization of energy (^|) over u k . Substi- 
tuting these u£ into Eq. (|5|) we obtain the effective Hamiltonian H = H c + Hd where the deformational interaction 
Hd can be written as 

H * = S ( n - kf ik + E "VST J G k U k f lk + ff kf 2 R kj = ^2 + H d3 + H di . (36) 



6 



Here the matrix Gk = G k ^ is inverse to the force constant matrix A^f , and the matrix product aBc means the sum 
a-aB a pb/3. The term Hd2, Hd3 and Hdi in (36) describes the pairwise, three-particle and four-particle deformational 
interaction, respectively: 



H, 



,12 



2 ^ 



1 



n(r)$ 2 (r - r')n(r'), = - ^ £ n(r)$ R (r - r>(r')n(r' + R), 



(37) 



R 



ff« = ^J n(r)n(r + R)$f' R '(r - i>(r')n(r' + R'), (38) 

r,r' ,R,R' 



where the potential <& 2 , or $ RR is given by the expression: 

{$ 2 (r) ; $ R (r); $ RR '(r)} e* kr {^ k G k fi k ; fi* k G k f R + f 2 ^G k f lk ; f R *G k f R k } . (39) 

k 

As the matrix G k in (39) at small k includes the well-known 'elastic singularity '0: G k ~ 1/fc 2 , each of terms Hd2, 
Hd3 and H44 in ( |36| ) includes the long-ranged elastic interaction. The formation of twinned structures discussed below 
is determined by the four-particle interaction Hdi- The rest deformational interactions, Hd2 and Hd3, for the single- 
phase Llo ordering under consideration lead just to some quantitative renormalizations of chemical interaction H c in 
(||) which are usually small and insignificant. Therefore, below we retain in the deformational interaction ( |36| ) only 
the last term Hdi- Let us also note that each term in the sum ( |38| ) for Hdi at fixed r and r' has the order of magnitude 
of the above-mentioned ordering-induced elastic energy e°l d which usually is small. Thus the interaction HdA can be 
significant only because of 'coherent' contributions of many sites r and r' due to the long-ranged elastic interaction. 
Therefore, local fluctuations of occupations.*/!?), in the interaction Hdi are insignificant, and it can be treated in the 
'kinetic mean-field approximation' (KMFA)t5~tJ which neglects such fluctuations and corresponds to the substitution 
in (|3^) of each occupation operator n(r) Va2-.its.mean value c(r) = (n(r)) where (. . .) means averaging over the space- 
and time-dependent distribution functiont^l~t£l. Therefore, in considerations of Al— >Llo transformations below we 
approximate the total effective Hamiltonian H in (j5|) by the following expression: 

H = H c + H di = H c {n(v)} + ^ ^ c(r)c(r + R)<I> R ' R '(r — r')c(r')c(r' + R') (40) 

r,r',R,R' 

where the potential $ R ' R (r) is given by the last equation ((3 



III. MODELS AND METHODS OF SIMULATION 

To-simuLate Al^Llo transformations in an alloy with the Hamiltonian ( ^p| ) we use the methods described in 
RcfsJi3 ancu3 to be referred to as-4, and II, respectively. Evolution of atomic distributions is described by the kinetic 
tetrahedron cluster field methodt-J in which mean occupations = c(rj) = (n(rj)) averaged over the space- and 
time-dependent distribution function obey the kinetic equation (1. 10): 

dc t /dt = 2J2 Mij sinh[/3(A i - Aj)/2]. (41) 
j 

Here (3 = 1/T is the inverse temperature; My is the generalized mobility proportional to the configurationally 
independent factor j nn in the probability of an inter-site atomic exchange Az «-> Bj between neighbouring sites i 
and j per unit time; and A^ =.A*{c,,} is the local chemical potential equal to the derivative of the generalized free 
energy F{c{\ defined in RefsJiMj with respect to \ = dF/dci. The expression for My = Myjcfc} employed in 
our simulations is given by Eq. (1.12) with the asymmetrical potential Uj taken zero for simplicity, while the local 
chemical potential A^ now is the sum of the chemical and the deformational term, A^ and Xf . The microscopical 
expressions for A^ are given by equations (1.13) - (1-16) which include only chemical interactions = v£j, while 
the deformational contribution Af = Xf (r^) is the variational derivative of the second term in ( fic| ) with respect to 
Cj = c(rj): 

Xf (r) = SH di /Sc(r) = 2 ^ c(r + R)$f R ' (r - r')c(r')c(r' + R'). (42) 

r'.R.R' 
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For the chemical interaction we employ the five alloy models used in I and II: 

1. The second- neighbour interaction model with the nearest-neighbour interaction V\ — 1000 K (in the Boltzmann 
constant fcs units) and V2/V1 = e = —0.125. 

2. The same model with e = —0.25. 

3. The same model with e = —0.5. r— ■ 

4. The fourth-neighbour interaction model with v n estimated by Chassagne et alE3 from their experimental data 
for disordered Ni-Al alloys: v 1 = 1680 K, v 2 = -210 K, v 3 = 35 K, and v 4 = -207 K. 

5. The fourth-neighbour interaction model with v\ — 1000K, v-ijv\ — —0.5, Vs/v% = 0.25, and v^/vi = —0.125. 

The effective interaction range Ri n t for these models monotonously increases with the model number. Therefore, a 
comparison of the simulation results for these models enables one to study the influence of Rint on the microstructural 
evolution. The critical temperature T c for the phase transition Al— >Llo in the absence of deformational interaction 
HdA (which seems to-have little effect on T c in our simulations) for model 1, 2, 3, 4 and 5 is 614, 840, 1290, 1950 and 
2280 K, respectively^. 

For the Kanzaki force fv^. entering the expression ([39|) for the potential <J>^' R (r) in fl42| ) we use Eqs. ( |34| ) and 
(jsBj). The interaction parameters ip a and in these equations can be expressed via spontaneous deformations u 
and e using Eqs. ( p9f ) and (|3^). For simplicity we assume the volume striction to be small (as it is for the Co-Pt 
alloys mentioned above): u ~ 0, while the tetragonal distortion will be characterized by its maximum value e m in a 
stoichiometric alloy, i. e. by the value e in (|2^) at r\\ = 0.5. Therefore, interactions ip a and (pb in our simulations are 
determined by the relations: 

tp a = -a(c n - ci 2 )e m /3; tp b = a(c n - c 12 )s m /12. (43) 

For the lattice constant a in (^) we take a-typical value a ~ 4 A, and for the elastic constant (en — C12), the value 
0.97 Mbar corresponding to FCC platinumEj. 

For the force constant matrix Ak (which determines the matrix Gk = (Ak) -1 in Eq. (|39"|)) we use the model 
described in Refs.Ej'Ej. It corresponds to a Born- von Karman model with the first- and second-neighbour force con- 
stants only, and the second-neighbour constants are supposed to correspond to a spherically symmetrical interaction. 
This model includes three independent force constants whisk are expressed in terms of elastic constants Cik , and these 
constants were chosen equal to those of the FCC platinumEa: cu = 3.47 Mbar, C12 = 2.5 Mbar, and C44 = 0.77 Mbar. 

As it was discussed in I and II, the transient partially ordered alloy states can be described using either mean 
occupations c, = c(r^) or local order parameters r}^ and local concentrations Cj defined by Eqs. (1.24) and (1.25). The 
simulation results below are usually presented as the distributions of quantities 77? = 77^ + 77^ + 77^ , to be called the 
'77 -representation', and these distributions are similar to those observed in the experimental transmission electron 
microscopy (TEM) imagesEj. 

Our simulations were performed in the FCC simulation boxes of sizes Vj, — L 2 x H (where L and H are given 
in the lattice constant a units) with periodic boundary conditions. We used both quasi-2D simulations with H = 1 
and 3D simulation with H ~ L. For the given coordinate z = na (with n = for 2D simulation) each of figures 
below shows all FCC lattice sites lying in two adjacent planes, z = na and z = (n + l/2)a. The point (x,y) with 
(x/a,y/a) equal to (l,m), (I + 1/2, m), (I + 1/2, m+ 1/2) or (l,m + 1/2) in the figures corresponds to the lattice site 
with (x/a, y/a,z/a) equal to (l,m,n), (I + 1/2, m, n+ 1/2), (I + 1/2, m + 1/2, n) or (l,m+ 1/2, n + 1/2), respectively. 
Therefore, at V& = L 2 x H the figure shows AL 2 lattice sites. 

The simulation methods were the same as in I and II. In simulations of Al^Llo transformation the initial as- 
quenched distribution Cj(0) was characterized by its mean value c and small random fluctuations Scf, usually we used 
Sci = ±0.01. The distribution of initial fluctuations Sci for the given simulation box volume VJ, was identical for all 
models and the same as that used in II. The sensitivity of simulation results to variations of these initial fluctuations 
Sci was discussed in II and was found to be insignificant for the features of evolution discussed below. 

IV. KINETICS OF A1^L1 TRANSFORMATIONS IN SYSTEMS WITH AN EXTENDED OR 

INTERMEDIATE INTERACTION RANGE 

As discussed in I, II and below, features of microstructural evolution under Al— >Llo and Al^Ll2 transitions 
sharply depend on the effective interaction range Rmt in an alloy. In this section we discuss Al— >Llo transitions for 
the systems with an extended or intermediate interaction range, such as our models 5 and 4, while the short-range- 
interaction systems are considered in the next section. 
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Some results of our simulations are presented in figures 1-8. The symbol A or A in these figures corresponds to an 
Llo-ordered domain with the tetragonal axis c along (100) and the positive or negative value, respectively, of the order 
parameter 771 ; the symbol B or B, to that for the c-axis along (010) and the order parameter r]2] and the symbol C or 
C, to that for the c-axis along (001) and the order parameter 773. Figure || shown in the c-representation illustrates 
the occupation of lattice sites for each domain type. The APB separating two APDs with the same tetragonal axis 
(i. e. APDs A and A, B and B or C and C) will be for brevity called the 'shift-APB', and the APB separating the 
APDs with perpendicular tetragonal axes will be called the 'flip- APB'. 

Before discussing figures 1-8 we remind the general ideas about the formation of twinned structures [1-11]. To 
avoid discussing the problems of nucleation, in this work we consider the transformation temperatures T lower 
the ordering spinodal temperature T s . Then the evolution under Al— >L1 transition includes the following stages 

(i) The initial stage of the formation of finest Llo-ordered domains when their tetragonal distortion makes still little 
effect on the evolution and all six types of APD are. present in microstructures in the same proportion. It corresponds 
to the so-called 'mottled' contrast in TEM imagesBfl 

(ii) The next, intermediate stage which corresponds to the so-called 'tweed' contrast in TEM images. The tetragonal 
deformation of the Llo-ordered APDs here leads to the predominance of the (llO)-type orientations of flip-APBs, but 
all six types of APD (i. e. AP-Ds with all three orientations of the tetragonal axis c) are still present in microstructures 
in comparable proportionsou. 

(iii) The final, polytwinned stage when the tetragonal distortion of the Llo-ordered APDs becomes the main factor 
of evolution and leads to the formation of (llO)-type oriented twin bands. Each band includes only two types of APD 
with the same c axis, and these axes in the adjacent bands-are 'twin' related, i. e. have the alternate (100) and (010) 
orientations for the given set of the (llO)-oriented bandsoQ. 

The thermodynamic driving force for the (llO)-type orientation of flip-APBs is the gain in the elastic energy 
of adjacent APDs: at other orientations this energy increases under the growth of an APD proportionally to its 
volumeETHl. For an APD with the characteristic size I and the surface Sd, this elastic energy E v el ~ cJe 2 Sdl begins 
to affect the microstructural evolution when it becomes comparable with the surface energy E s ~ aSd where a is the 
APB surface tension. The 'tweed' stage (ii) corresponds to the relation E v el ~ E s or to the characteristic APD size 

l ~ a/c_e 2 , (44) 

and so this size sharply increases under decreasing distortion e. 

Figures 1-7 illustrate quasi-2D simulations for which microstructures include only edge-on APBs normal to the 
(001) plane. The elimination of the volume-dependent elastic energy mentioned above is here possible only for the 
(100) and (OlO)-oriented APDs separated by the (110) or (llO)-oriented APBs, while in the (OOl)-oriented APDs 
C and C this elastic energy is always present. Therefore, the tweed stage (ii) in these simulations corresponds to 
both the predominance of (110) or (llO)-oriented APBs separating domains A or A from B or B and the decrease 
of the portion of domains C and C in the microstructures. In the 3D case each of three posible types of a polytwin, 
that without (001), (100), or (OlO)-oriented APDs, can be formed in the given part of an alloy stochastically due to 
the local fluctuations of composition [1-7]. It is illustrated, in particular, by 3D simulation shown in figure |8j, while 
quasi-2D simulations describe the formation of only one polytwin type mentioned above. 

The distortion parameter |e m | = 0.1 for the simulations shown in figures |l|-^ was chosen so that the APD size 
lo in Eq. (H) characteristic for manifestations of elastic effects has the scale typical for real CoPt-type alloys. In 
particular, if we take a conventional assumption that the APB energy a is proportional to the transition temperature 
T c : ~ T c f(T') where / is some function of the reduced temperature T' = T/T c , then using the relation e m = e/A-qf 
and the parameters e, 771 , and T c for CoPt and for our models mentioned above we find that the right-hand side of Eq. 
(Eij) for models 5 and 4 at |e m | = 0.1 is close to that for the CoPt alloy at similar T" values within about ten percent. 
Therefore, the microstructures at both the initial stage (i) and the tweed stage (ii) can be reproduced by figures [l]-|| 
with no significant distortion of scales. Under a furher growth of an APD its size I becomes comparable with the 
simulation box size L, and the periodic boundary conditions begin to significantly affect the evolution. Therefore, the 
later stages of transformation can be more adequately simulated if we reduce the characteristic size Iq in Eq. (|44) 
using the larger values of theparameter e m , such as \e m \ = 0.15 — 0.2 used in the simulations shown in figures |f^f 

Let us first discuss figures corresponding to a 'realistic' distortion parameter \e m \ — 0.1. The initial stage (i) 
in these figures corresponds to frames |l](a)-jl](b), ||(a)-0(b), and 0(a)-^(b)- the tweed stage (ii), to frames |l|(c)-|l|(d), 
||(c)-||(e), and ||(c); and the twin stage (iii), to frames |l](e)-|l](f),||(f), and||(d). |— , 

The detailed consideration of the initial stage for models 4 and 5 neglecting the deformational effectal3 revealed 
the following features of evolution: 

(a) The presence of abundant processes of fusion of in-phase domains which are one of main mechanisms of domain 
growth at this stage. 

(b) The presence of peculiar long-living configurations, the quadruple junctions of APDs (4-junctions) of the type 
A1A2A1A3 where A2 and A3 can correspond to any two of four types of APD different from Ai and Ai. 







(c) The presence of many processes of 'splitting' of a shift-APB into two flip- APBs which lead to either a fusion 
of in-phase domains mentioned in point (a) (s — > / process) or a formation of a 4-junction mentioned in point (b) 
(s — > 4j process). 

Figures [l]-0 show that all these microstructural features are also present when the deformational effects are taken 
into account, and not only at the initial stage (i) but also at the tweed stage (ii). In particular, the bcginningand the 
end of an s — ► / process (marked by the single and the thick arrow, respectively) can be followed in frames 0(a) and 
0(b); 0(c) and 0(d) ; 0(d) and 0e; 0(b) and |(c); §(c) and |(d); and |(a) and |(b). ThcJusion with the disappearance 
of an intermediate APD which initially separates two in-phase domains to be fusedtj can be followed in the lower 
right part of frames [l](a) and 0(b) and in the upper right part of frames 0(b) and 0(c) (which is marked by a thick 
arrow in frames 0(b) and0(c), respectively). A number of long- living 4-junctions marked by thin arrows are seen in 
frames 0(a) -0(d) , |2|(a)-g(c), and 0(a). An s — ► 4j process can be followed in the lower right part of frames 0(a)-0(c) 
. The processes and configurations (a), (b) and (c) can also be seen in figures 4-7 discussed below—. 

Frames 0(a)-0(e) also display some (lOO)-oriented and thin conservative APBs. As discussed iiiEil and below, such 
APBs are most typical of the short-range-interaction systems where they have a low surface energy (being zero for 
the stoichiometric nearest-neighbor interaction model) unlike other, non-conservative APBs. Under an increase of the 
interaction range, as well as temperature or the deviation from stoichiometric composition, the anisotropy in the APB 
surface energy sharply decreasesE3. Therefore, in figure (and figure [5] below) corresponding to the intermediate- 
range-interaction model 4 the conservative APBs are few but observable, while for the extended-range-interaction 
model 5 in figure 0, as well as at elevated T or significant 'non-stoichiometry' 5c = (0.5 — c) in figures or for model 
4, such APBs are absent entirely. 

Comparison of figures and illustrates the sharp dependence of microstructural evolution on the transformation 
temperature T. Under elevating this temperature to values near the critical temperature T c : (T c — T) < 0.1 T c , both 
flip and shift-APBs notably thicken, the anisotropy in their surface energy falls off, while the characteristic size of 
initial APDs (formed after a rapid quench Al— >Llo) increases. The latter is related to an increase at T — > T c of 
the characteristic wavelength for the ordering instability which is due to the narrowing of the interval of effective 
wavenumbers q = k — k s near the superstructure vector k s for which the ordering concentration waves are unstable 
at T < T c . 

Frames 0(c)-0(d), 0(c)-0(e), and 0(c) show evolution at the tweed stage. They illustrate, in particular, kinetics 
of the (llO)-type alignment of APBs between APDs A or A and B or B at this stage, as well as a "dying out" of 
(lOO)-oriented APDs C and C. These frames also show that in the simulation with a realistic distortion parameter 
\e rn \ = 0.1 (fitted to the structural data for CoPt) the APD size l Q) characteristic of the tweed stage is About 
(20 — 40) a. It agrees with the order of magnitude of this size observed in the CoPt-type alloys FePt and FePdEl u. 

As mentioned, the final, twin stage of the transformation can be more adequately simulated with the larger values 
of parameter \e m \ which are employed in the simulations shown in figures 0-0. Before discussing these [figures we 
note some typical configurations observed, in experimental studies of transient twinned microstructuresoQ seen, for 
example, in figures 5, 6, 9 and 2 in Rcfsffl,Q andlia, respectively: 

(1) semi-loop-like shift-APBs adjacent to the twin band boundaries; 

(2) 'S-shaped' shift-APBs stretching across the twin band; and 

(3) short and narrow twin bands (for brevity to be called 'microtwins') which lie within the larger twin bands and 
usually have one or two shift-APBs near their edges. 

Comparing our results with experiments one should consider that due to the limited size of the simulation box 
the twin band width d in our simulations has the same order of magnitude as the APD size Iq (|44|) characteristic 
of the tweed stage, while in experiments d usually much exceeds loan. Therefore, the distribution of shift-APBs 
within twin bands in our simulation is usually much more close to equilibrium than in experiments. In spite of this 
difference, the simulations reproduce all characteristic transient configurations (l)-(3) and elucidate the mechanisms 
of their formation. In particular, both the semi-loop and the S-shaped shift-APBs are formed from regular-shaped 
approximately quadrangular APDs (characteristic of the beginning of the twin stage) due to the disappearance 
of adjacent APDs which are 'wrongly-oriented' with respect to the given twin band. The formation of semi-loop 
configurations is illustrated by frames 0(d)- 1 f) , B|(d)-||(e) , and 7(b)-0(c); while the formation of S-shaped APBs 
can be seen in frames 0(d)-[ff), l(d)-lff), ||d)-|(e), and 0(b)-f|(c). The formation and evolution of microtwins 
is illustrated by frames 0(c)- |(d) and ||(c)- 5(d). These frames show that the microtwin is actually a small and 
narrow APD for which deformational effects are strong enough to align its flip- APBs along the (llO)-type directions. 
However, the standard mechanism of coarsening via the growth of larger APDs at the expence of smaller ones leads 
to the shrinking and eventually to the disappearance of a microtwin which is usually accompanied by the formation 
of S-shaped and/or semi- loop shift-APBs. The latter is illustrated by frames 0(d)-|4(f) and 0(d)-0(e). Let us also 
note that the microtwin configuration shown in frame 0(d) is strikingly similar to that seen in the central part of 
experimental figure 2 in Ref.EJ. 

Let us now discuss the final, 'nearly equilibrium' microstructures shown in last frames of figures 0-0 A characteristic 
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feature of these microstructures is a peculiar alignment of shift- APBs: within a (lOO)-oriented twin band in a (110)- 
type polytwin the APBs tend to align normally to some direction n = (cos a, sin a, 0) characterized by a 'tilting' angle 
a between the band orientation and the APB plane. Figures [l|-|| show that this tilting angle is not very sensitive 
to the variations of temperature or concentration but it sharply depends on the interaction type, particularly on the 
interaction range. For the extended-range- interaction model 5 this angle is close to 7r/4 (slightly exceeding this value), 
while for the intermediate-range- interaction model 4 it is notably less than tt/4. A similar alignment of APBs for the 
short-range interaction systems is illustrated by figures p|-|Tl| where the tilting angle is close to zero. 

A phenomenological theory of this interaction-dependent tilting of APBs within nearly-equilibrium twin bands has 
been suggested inlill. The tilting is explained by the competition between the anisotropy of the APB surface tension 
a and a tendency to minimize the total APB area within the given twin band which corresponds to a = tt/4. For 
the alkw-systems with both the intermediate and the short interaction range the surface tension a (a) is minimal at 
a = OlililS. Thus for such alloy systems the tilting angle is less than 7r/4, and it decreases with the decrease of-thc 
interaction range. For the extended-range-interaction systems the anisotropy of the APB surface tension is weakE§E3, 
and so the .tilting angle is close to 7r/4. Therefore, the comparison of experimental tilting angles with theoretical 
calculationsO can provide both qualitative and quantitative information on the effective chemical interactions in an 
alloy. 

The alignment of shift- APB|S. discussed above seems to be clearly seen in the experimental microstructure for 
CoPt shown in figure 5 of Refill where the tilting angle is notably less than 7r/4. It can indicate that the effective 
interactions in CoPt have an intermediate interaction range. This agrees with the usual estimates of these interactions 
for transition metal alloys, see e.g.li3c3'E3. 

Comparison of figures and u\ illustrates the influence of temperature T on the evolution. Under elevating T 
we again observe a thickening of APBs, as well as a coarsening of initial APDs. Frames 0(d)-|7](f) also illustrate a 
process of 'transverse coarsening' of twin bands via a shrinkage and disappearance of some microtwinned bands. Such 
transversepGoarsening appears to be seen in a number of experimental microstructures, for example, in figure 6 intl or 
figure 2 int3. Frames 0(d)-Q(f) show that the thermodynamic driving force for such transverse coarsening is mainly 
the gain in the surface energy of shift- APBs due the decrease of their total area under this process. 

Figures || and || illustrate the concentration dependence of the evolution. The non-stoichiometry 5c — (0.5 — c) 
affects the evolution similarly to temperature T: under an increase of both 5c and T all APBs thicken, while shift- 
APBs become less stable with respect to flip-APBsE3. The latter leads to an enhancement of processes of splitting of 
shift- APBs as well as of the transverse coarsening mentioned above; it is illustrated by frames ^(b)- |](f). 

Some results of a 3D simulation with V& = 52 2 x 30 are presented in figure ||. In this figure we employ the o 
representation (described in the caption) in which the regions containing the vertical or horizontal lines (that is, 
the vertical or horizontal crystal planes filled by A atoms) correspond to the APDs with the (100) or (OlO)-oriented 
tetragonal axis, respectively, while the checkered regions correspond to the APDs with the tetragonal axis normal to 
the plane of figure. This simulation aimed mainly to complement 2D simulations with an illustration of geometrical 
features of 3D microstructures. Figure || illustrates, in particular, a stochastic formation of different polytwin sets 
with three possible types of orientation mentioned above. A limited size of the simulation box prevent us from a 
detailed consideration of evolution with this 3D simulation. Therefore, below we discuss only the problem of a 3D 
orientation of tilted shift- APBs in final, 'nearly-equilibrium' microstructures. 

Let us consider a (lOO)-oriented twin band in the form of a plate of height h, length I, and width d in the direction 
(001), (HOL and (110), respectively, with d < h -C I (which is a typical form of twin bands observed in TEM 
experimentsoiH) . The equilibrium orientation of a plane shift- APB in this band corresponds to the minimum of its 
energy E s = aS where S is the APB area and a is the surface tension determined mainly Jsy the angle a between the 
APB orientation n = (sin a, cos a cos ip, cos a sin ip) and the band orientation no = (100)til. For the 'needle-shaped' 
twin band under consideration the upper and the lower boundary of a shift-APB usually lies at the top and the 
bottom edge of this band, respectively. Minimization of energy E s in this case yields: ip = 0, i. e. the APB is normal 
to the (0011-plane, and its orientation n = (sin a, cos a, 0) is determined by the interaction-dependent tilting angle 
a defined intil. This conclusion seems to be supported by the present 3D simulation: the lower and the upper tilted 
shift-APB within the (OlO)-oriented twin band below the main diagonal of frame @(c) corresponds to the grey line 
stretching across the checkered region in frame ||(e) and |](f), respectively, and both these lines are approximately 
normal to the (001) plane. 

V. KINETIC FEATURES OF A1->L1 TRANSFORMATIONS IN THE SHORT-RANGE INTERACTION 

SYSTEMS 

As mentioned, transient microstructures under Llo ordering for the short-range interaction systems include many 
conservative APBs. Such APBs are virtually immobile, and so the evolution is realized via motion of other, non- 
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conservative APBs which results in a number of peculiar kinetic featurestJ'ta. The-initial stage of the Al^Llo 
transformation for the short-range interaction systems was discussed in detail in Ref.ES. In this section we consider 
the tweed and twin stages of such transformations and note the differences with the case of systems with the larger 
interaction range. 

Some results of our simulations for the short-range-interaction systems are presented in figures 10 In these 
simulations we used sufficiently high temperatures T" > 0.9 — 0.8 to accelerate evolution to final, 'nearly-equilibrium' 
configurations as the presence of conservative APBs slowes down (or even 'freezes') this evolution, particularly at low 
X". 

Figure || illustrates the evolution for model 1; as discussed inEl, this model seems to correspond to the Cu-Au-type 
alloys. A distinctive feature of microstructures shown in figure]^ is a predominance of the above-mentioned conservative 
APBs with the (lOO)-type orientation. Frames ||(a)-^(c) show both the conservative shift- APBs (c-shift APBs) and the 
conservative flip- APBs (c-flip APBs) also illustrating their orientational properties!^. For quasi-2D microstructures 
with edge-on APBs shown in figure |9|, c-shift APBs separating APDs A and A (c-APBs A-A) are horizontal; c-APBs 
B-B are vertical; and c-APBs C-C can be both horizontal and vertical; c-flip APBs (A or A)-(C or C) (which separate 
APDs A or A from C or C) are horizontal; c-APBs (B or B)-(C or C) are vertical; and c-APBs (A or A)-(B or B) 
should lie in the plane of figure and thus they are not seen in figure ^. Figure ^ also shows that the conservative APBs 
are notably thinner than non-conservative ones, particularly so for c-flip APBs. 

Frames S(a)-p|(c) show that at first stages of evolution the portion of conservative APBs with respect to non- 
conservative ones increases, due to the lower surface energy of the c-APBs. Later on, with the beginning of the tweed 
stage, the deformational effects become important leading to a dying out of both APDs C and C and their c-flip APBs. 
However, the conservative shift- APBs within twin bands survive, and in the final frame ^d they are mostly 'step-like' 
consisting of (lOO)-type oriented conservative segments and small non-conservative ledges. These step-like APBs can 
be viewed as a 'facetted' version of tilted APBs discussed above and seen in figures [l|-||. Such step-like APBs were 
observed under the Llo ordering of CuAu aed some CuAu-based alloysQ, and they pare also similar to those observed 
under the LI2 ordering in both simulations!^ and experiments for the CU3AU alloycJ. 

As it was repeatedly noted in Ref.Ej and above, an increase of non-stoichiometry 5c — (0.5 — c) or temperature 
T leads to a sharp decrease of both the anisotropy of the APB energy and the energy preference of conservative 
APBs with respect to non-conservative ones. Therefore, under an increase of 6c or T the portion of conservative 
APBs in transient microstructures falls off, and at sufficiently high 6c or T such APBs are not formed under the 
transformation at all. It results in drastic microstructural changes of evolution, including sharp, phase-transition-like 
changes in morphology of aligned shift- APBs within twin bands, from the 'faceting' to the 'tilting'. This is illustrated 
by figure |l0| which shows the evolution of model 1 at a significant non-stoichiometry 6c = 0.06, and this evolution is 
qualitatively different with that for a stoichiometric alloy shown in figure ^. 

Figure illustrates the transition from the 'facetted' to the 'tilted' morphology of shift-APBs within nearly- 
equilibrium twin bands under variations of T or 6c for model 2. An examination of intermediate stages of transfor- 
mations illustrated by this figure shows that the morphological changes are realized via some local bends of facetted 
APBs. It is also illustrated by a comparison of frames 0(a), 0(c) and 0(d) with each other. Therefore, the 'mor- 
phological phase transition' mentioned above is actually smeared over some interval of temperature or concentration. 
However, frames 0(a)-0(d) show that the 'intervals of smearing' of such transitions can be relatively narrow. 



VI. CONCLUSION 

Let us summarize the main results of this work. The earlier-described master equation approacbEl'O is used to 
study the microstructural evolution under Llo-type orderings in alloys, including the formation of twinned structures 
due to the spontaneous tetragonal deformation inherent to such orderings. To this end we first derive a microscopical 
model for the effective interatomic deformational interaction which arise due to the so-called Kanzaki forces describing 
interaction of lattice deformations with site occupations. This model generalizes an analogous model of Khachaturyan 
for dilute alloyal3 to the physically interesting case of concentrated alloys. We take into account the non-pairwise 
contribution to Kanzaki forces, and the resulting effective interaction Hd is non-pairwise, too, unlike the case of 
dilute alloys. This effective interaction describes, in particular, the lattice symmetry change effects under phase 
transformations, such as the tetragonal distortion mentioned above. Assuming the non-pairwise Kanzaki forces to be 
short-ranged, we can express the deformational interaction Hd in terms of two microscopical parameters which can 
be estimated from the experimental data about the lattice distortions under phase transformations. We present these 
estimates for alloys Co-Pt for which such structural-data are availableEd. 

Then we employ the kinetic cluster field methodEjIlS to simulate Al— >Llo transformation after a quench of an alloy 
from the disordered Al phase to the single-phase Llo field of the phase diagram in the presence of deformational inter- 
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action Hd- We consider five alloy models with different types of chemical interaction, from the short-range-interaction 
model 1 to the extended-range-interaction model 5, at different temperatures T, concentrations c, and spontaneous 
tetragonal distortions e. We use both 2D and 3D simulations, and all significant features of microstructural evolution 
in both types of simulation were found to be similar. 

The evolution under Al— >Llo transition can be divided into three stages, in accordance with an increasing im- 
portance of the defocmational interaction Hd'- the 'initial', 'tweed' and 'twin' stage. For the initial stage (discussed 
in detail previouslyE3) , the deformational effects are insignificant. For the tweed stage, the effects of Hd become 
comparable with those of chemical interaction H c and lead to the formation of specific microstructures discussed in 
section 4. For the final, twin stage the tetragonal distortion of Llo-ordered antiphase domains (APDs) becomes the 
main factor of the evolution and leads to the formation of (llO)-type oriented twin bands. Each band includes only 
two types of APD with the same tetragonal axis, and these axes in the adjacent bands are 'twin' related, i. e. have 
the alternate (100) and (010) orientations for the given set of (llO)-type oriented bands. 

The microstructural evolution strongly depends on the interaction type, particularly on the interaction range Ri n t- 
For the systems with an extended or intermediate Rint at both the initial and the tweed stage we observe the following 
features (mentioned previousljtj for the initial stage): (a) abundant processes of fusion of in-phase domains; (b) a 
great number of peculiar long-living configurations, the quadruple junctions of APDs described in section 4; and 
(c) numerous processes of 'splitting' of an antiphase boundary separating the APDs with the same tetragonal axis 
('shift- APB') into two APBs separating the APDs with perpendicular tetragonal axis (flip-APBs). The simulations 
also illustrate a sharp temperature dependence of the evolution, in particular, a notable increase of both the width 
of APBs and the characteristic size of initial APDs under elevating T. The deviation from stoichiometry affects the 
evolution similarly to temperature: under an increase of both non-stoichiometry 5c = (0.5 — c) and T all APBs thicken, 
while shift-APBs become less stable with respect to flip-APBs. 

For the twin stage, our simulations reveal the following typical features of transient microstructures: (1) semi-loop- 
like shift-APBs adjacent to the twin band boundaries; (2) 'S-shaped' shift-APBs stretching across the twin band; 
(3) short and narrow twin bands ('microtwins') lying within the larger twin bands; and (4) processes of 'transverse 
coarsening' of twinned structures via a shrinkage and disappearance of some microtwins. All these features agree with 
experimental observationsoQ. For the final, ncarly-cquilibrium twin bands the simulations demonstrate a peculiar 
alignment of shift-APBs with a certain tilting angle between the band orientation and the APB plane, and this 
tilting angle sharply depends on the interactiomJype, particularly on the interaction range Rint- Such alignment 
of APBs seems to be observed in the CoPt alloyEl, and a comparison of experimental tilting angles with theoretical 
calculations^ can provide information about the effective interactions in an alloy. 

A distinctive feature of evolution for the short-range-interaction systems is the presence of many conservative 
APBs with the (100)-type orientation. The conservative flip-APBs disappear in the course of the evolution, but 
the conservative shift-APBs survive and are present in the final twinned microstructures. Such 'nearly equilibrium' 
shift-APBs are mostly 'step-like' consisting of (100)-type oriented conservative segments and small non-conservative 
ledges, which can be viewed as a 'facetted' version of tilted APBs mentioned above. This (100)-type alignment of 
shift-APBs within twin bands seems to agree with ava ilab le experimental observations for the CuAu alloyQ for which 
chemical interactions are supposed to be short-rangedc3ll3. 

Under an increase of non-stoichiometry 5c or temperature T the energy preference of conservative APBs with respect 
to non-conservative ones decreases, and the portion of conservative APBs in the microstructures falls off. It results in 
drastic microstructural changes, including sharp, phase-transition-like changes in morphology of aligned shift-APBs 
within twin bands, from their 'faceting' to the 'tilting'. Such 'morphological phase transitions' are actually smeared 
over some intervals of temperature or concentration, but the simulations show that the intervals of smearing can be 
narrow. 

Finally, let us make a general-remark about kinetics of multivariant orderings in alloys, such as the LI2, Llo and 
DO3 orderings discussed in RefsJl3'E3E3 and in this work. It is known that the thermodynamic behavior of different 
systems under various phase transitions reveals features of universality and insensitivity to the microscopical details 
of structure, particularly in the critical region near thermodynamic instability points. The results of this and other 
studies of multivariant orderings show that such universality does not seem to hold for their phase transformation 
kinetics, at least outside the critical region (which for such orderings is usually either quite narrow or absent at all). 
The microstructural evolution reveals a great variety of peculiar features, the detailed form of which sharply depends 
on the type of interatomic interaction, the type of the crystal structure and ordering, the degree of non-stoichiometry, 
and other 'non-universal' characteristics. 
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FIG. 1. Temporal evolution of the extended-interaction-range model 5 under the phase transformation Al^Llo shown in 
the ^-representation for the simulation box size V b = 128 2 x 1 at the maximum tetragonal distortion parameter |e m | = 0.1, 
c = 0.5, the reduced temperature T' — T/T c = 0.7, and the following values of the reduced time t' = tj nn : (a) 10; (b) 20; 
(c) 50; (d) 100; (e) 250; and (f) 280. The grey level linearly varies with r/f — 77^ +ri2i + 7 ?fi between its minimum and maximum 
values from completely dark to completely bright. The symbol A, A, B, B, C or C indicates the type of the ordered domain as 
described in the text. The thick, the thin and the single arrow indicates the fusion-of-domain process, the quadruple junction 
of APDs, and the splitting APB process, respectively, discussed in the text. 




FIG. 2. As figure [l| but for the intermediate-interaction-range model 4 at T' = 0.67 and the following values of t': (a) 10; 
(b) 20; (c) 50; (d) 100; (e) 250; and (f) 500. 



15 



FIG. 3. As figure |, but at T 1 = 0.92 and the following values of t': (a) 10; (b) 50; (c) 100; and (d) 200. 





FIG. 5. As figure |, but at |e m | = 0.15 and the following values of t': (a) 10; (b) 20; (c) 50; (d) 170; (e) 200; and (f) 300. 



FIG. 6. As figure |, but at c = 0.44 and the following values of t': (a) 10; (b) 20; (c) 50; (d) 400; (e) 750; and (f) 1100. 
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FIG. 7. As figure but at |e m | = 0.2, T' = 0.88 and the following values of t': (a) 10; (b) 20; (c) 30; (d) 50; (e) 300; and 
(f) 400. 




FIG. 8. As figure but for 3D simulation with Vb = 52 2 x 30 at |e m | = 0.2 and shown in the 'c-representation': the grey 
level linearly varies with d between its minimum and maximum values from completely dark to completely bright. The three 
upper frames correspond to the plane z = 10 a and the following values of t': (a) 10; (b) 20; and (c) 325. The three lower 
frames correspond to t' = 325 and the following planes: (d) y = 0; (e) y = 10 a; and (f) y = 36 a. 
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FIG. 9. As figure [l] but for the short-range-interaction model 1 at |e m | = 0.15, c = 0.5, T' = 0.9 and the following values of 
t': (a) 30; (b) 40; (c) 60; and (d) 120. 



FIG. 10. As figure | but at c = 0.44 and the following values of t': (a) 110; (b) 140; (c) 200; and (d) 350. 
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